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Abstract 

The T2K long-basclinc neutrino oscillation experiment in Japan needs precise 
predictions of the initial neutrino flux. The highest precision can be reached 
based on detailed measurements of hadron emission from the same target as 
used by T2K exposed to a proton beam of the same kinetic energy of 30 GeV. 
The corresponding data were recorded in 2007-2010 by the NA61/SHINE 
experiment at the CERN SPS using a replica of the T2K graphite target. In 
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this paper details of the experiment, data taking, data analysis method and 
results from the 2007 pilot run are presented. Furthermore, the application 
of the NA61/SHINE measurements to the predictions of the T2K initial 
neutrino flux is described and discussed. 

Keywords: hadron production, long target, neutrino flux predictions 



1. Introduction 

Neutrino beams have become the major tool to perform studies of neu- 
trino properties. At the T2K long-baseline neutrino oscillation experiment 
in Japan [H 12] , a high-intensity neutrino beam is produced at J- PARC by a 
30 GeV proton beam impinging on a 90 cm long graphite target. Positively 
charged hadrons exiting the target (mainly vr and K mesons) are focused by 
a set of three magnetic horns and decay along a 96 m long decay tunnel. The 
flavour content and energy spectrum of the beam are measured at the near 
detector complex located 280 m away from the target station, and by the 
Super-Kamiokande (SK) detector at a distance of 295 km. For the first time 
in the history of accelerator-based neutrino experiments, T2K adopted the 
off-axis technique |3] to generate a dedicated neutrino beam with the off-axis 
angle set to 2.5° for both the near and far detectors. 

T2K was the first experiment to make a direct measurement of a non-zero 
value of the 6*13 mixing angle via — )■ z/g appearance. The published 90 % 
CL inclusion interval of 0.03(0.04) < sin^ 2^13 < 0.28(0.34) for the normal 
(inverted) mass hierarchy, 6cp = 0, sin^ 2^32 = 1 and Am^g = 2.4 x 10~^ eV^ 
was obtained with only 2 % of the final statistics P] . Later, these results were 
confirmed with greater precision by measurements of electron anti-neutrino 
disappearance at reactors t5j,iS]. With the same set of data T2K also provided 
new measurements of the neutrino oscillation parameters sin^ 2^32 and Am32 
by studying z/^ disappearance [7j , and aims at a precision of 1 % for sin^ 26'32 
and 3 % for 

Although neutrino beams provide well-defined and controlled sources of 
neutrinos, intrinsic uncertainties on the fluxes predicted with Monte Carlo 
(MC) simulations arise from models employed to simulate hadron emission 
from long nuclear targets used in accelerator-based experiments. In these 
typee of experiments, a non-negligible fraction of the neutrino flux actually 
arises from tertiary particles that are produced in hadronic re-interactions in 
the target. While a comparison of charged pion production from incident pro- 
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tons on cylindrical long and short targets has been published by the HARP 
Collaboration [8] , up to now, flux predictions have been (if ever) constrained 
by using dedicated hadron production measurements performed on thin nu- 
clear targets, e.g. HARP p+Al data t9j for K2K |10], HARP p+Be data [II] 
for MiniBooNE [[12] and SciBooNE, SPY p+Be data [13] for NOMAD [H]. 
T2K recently followed this approach by using the NA61/SHINE results on 
p+C interactions at 30 GeV extracted from measurements of hadron pro- 
duction in a thin (2 cm) graphite target [I5l HB] . 

Although such measurements provide constraints on the production of 
secondary particles in the primary interaction of the beam protons in the 
target, the lack of direct measurements of the production of tertiary parti- 
cles in re-interactions, and hence the use of sparse data sets to cover these 
contributions, limits the achievable precision of the flux prediction. The main 
motivation for measurements of hadron emission from a replica of the T2K 
target is therefore to reduce the systematic uncertainties on the prediction 
of the initial neutrino flux originating from products of interactions in the 
target. 

The NA61/SHINE (SPS Heavy Ion and Neutrino Experiment) exper- 
iment at the CERN Super Proton Synchrotron (SPS) is pursuing a rich 
physics program in various fields [ITJ [HI HU [20] from precise hadron pro- 
duction measurements for T2K and more reliable simulations of cosmic-ray 
air showers for the Pierre Auger and KASCADE experiments [211 [22], to the 
study of the properties of the onset of deconfinement with measurements of 
p+p, p+Pb and nucleus+nucleus collisions at the SPS energies. 

In addition to recently published thin-target (0.04 A/) measurements of 
charged pion and kaon production [T^ [TB] already used for the T2K neutrino 
flux predictions [U [7], the NA61 /SHINE collaboration studies hadron emis- 
sion from a replica of the T2K target (1.9 A/) exposed to a 30 GeV proton 
beam. A total of 0.2x10^ events were recorded during a pilot data taking 
in 2007. High statistics data were recorded in 2009 (4x10^ events) and 2010 
(10x10^ events). For the first time, the kinematical phase space of pions and 
kaons exiting the target and producing neutrinos in the direction of the near 
and far detectors is fully covered by a single hadron production experiment. 

The long-target analysis presented in this paper uses the low-statistics 
data collected in 2007. It however sets the ground for the ongoing analysis 
of high-statistics NA61/SHINE data with the replica of the T2K target. It 
demonstrates that high-quahty long-target data were successfully taken with 
the NA61/SHINE apparatus for T2K, and that such data can be used ef- 
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fectively to constrain the T2K neutrino flux predictions. A comparison of 
neutrino flux predictions based on thin-target hadron production measure- 
ments and long-target hadron emission measurements is performed as an 
illustration of the complete procedure. 

The paper is organised as follows: Section 2 briefly reviews the current 
T2K flux predictions based on the NA61/SH1NE thin-target data and points 
out the need for additional long-target measurements to improve the preci- 
sion of the predictions. Section 3 describes the NA61/SHINE experimental 
setup, kinematical coverage of the data, event selection and data normalisa- 
tion, reconstruction method and particle identification. The NA61/SHINE 
simulation chain is presented in Section 4. Yields of positively charged pions 
measured at the surface of the replica of the T2K target are given in Sec- 
tion 5. Possible strategies to use long-target measurements in the T2K beam 
simulation are proposed in Section 6 which also provides an illustration of 
the complete procedure. 

2. Requirements on hadron production data for the prediction of 
T2K neutrino fluxes 

The T2K beam MC simulation [2] is used to predict the initial neutrino 
flux at the near and far detectors. It comprises a full description of the beam 
line, including the target, magnetic horns, decay tunnel and beam dump. 
Hadronic interactions in the target are simulated by the FLUKA2008.3b [23] 
model. Measurements of particle emission from the replica of the T2K target 
are necessary to constrain the model calculations and to reach a 5 % precision 
on the absolute flux prediction as required by the T2K physics goals (i.e. 3 % 
precision on the ratio of the far to near fluxes for precision disappearance 
and Vf, appearance analyses). 

Predictions obtained for horn currents of 250 kA are shown in Fig. [l] for 
the z/^ and v^. fluxes at the near detector. The flux below 2 GeV predom- 
inantly (95 %) originates from the in-flight decay of positively charged pions 
focused by the magnetic horns of the beam line (see Ref. |2] for a detailed 
description of the T2K beam line). The flux is dominantly produced by 
the decay of positively charged kaons above 1.5 GeV, whereas at lower energy 
z/e's originate mostly from the decay of pions via the subsequent muon decay, 
i.e. 7r+ —7- [i^v^^ followed by /i"^ — e^Vey^l. Thus, pion production data can 
constrain most of the z/^ flux and a significant fraction of the v^, flux below 
2 GeV neutrino energy. 
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Figure 1: Prediction (based on FLUKA2008.3b and re-weighted by the NA61 thin target 
data) of the i/^ [left] and Vg [right] fluxes at the near detector of T2K. Different colours 
refer to the contributions of the various parent particles. 

In terms of hadron production measurements, neutrino fluxes can be de- 
composed into secondary and tertiary components. The secondary compo- 
nent originates from neutrino parents produced in the primary interaction 
of the beam protons in the target, e.g. secondary pions, p + C — )■ 7r+ + X. 
This secondary component can be constrained mainly by pion (and kaon) 
production cross-sections obtained from measurements on a thin target. The 
tertiary component refers to neutrino parents produced in interactions of sec- 
ondary particles, whether such interactions occur in the target or out of the 
target in the elements of the beam line. The contribution to the neutrino 
flux from parents produced in the target is therefore defined as the sum of 
the secondary component and the tertiary component due to interactions in 
the target. This contribution can be obtained from measurements of pion 
(and kaon) emission from a replica target. 

The dependence of the secondary and tertiary contributions on the neu- 
trino energy is depicted in Fig. |2] for the and Ve fluxes at the far detector. 
The secondary component contributes 60 % of the (z/g) flux at the peak 
of the beam energy spectrum (600 MeV). The remaining 40 % constitutes 
the tertiary component due to interactions in the target and elements of the 
beam line. Thus, thin-target measurements for T2K (i.e. positively charged 
pion and kaon inclusive production cross-sections at 30 GeV [Ul [I6j) can 
directly constrain up to 60 % of the i/^ (z/g) flux prediction. 

The lack of direct measurements of secondary interactions however re- 
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Figure 2: Secondary and tertiary components of the j/^ [top left] and [top right] 
fluxes at the far detector. The contribution of parents originating from the target sums 
up to 90 %, among which 60 % are due to the secondary component and 30 % due to 
re-interactions in the target (the in-target component). The relative contributions of the 
secondary and total in-target (secondary-|-tertiary in-target) components are shown for Vfj, 
[bottom left] and Ve [bottom right] as a function of energy. The dashed vertical line shows 
the location of the peak of the beam energy spectrum (600 MeV). Predictions are based 
on FLUKA2008.3b and re-weighted by the NA61 thin-target data. 



quires in most cases scaling to energies and nuclei relevant for the T2K ex- 
perimental setup, as well as extrapolating to uncovered regions of the kine- 
matical phase space. Such procedures have been used in addition for the 
T2K flux prediction. This brings in new sources of systematic uncertainties 
on top of the uncertainty of the measurements. 

As an example, the systematic errors of the z/^ and flux predictions at 
the far detector for the first published T2K analysis are depicted in Fig. |3] 
Details about the procedure developed to re-weight the original predictions of 
the T2K beam simulation (based on FLUKA2008.3b) with the NA61 thin- 
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target data are given elsewhere [21]. The total fractional error on the z/^ 
and z/e fluxes is about 15 % at the peak of the beam energy spectrum. At 
this energy the fractional error attributed to the re-weighting of tertiary 
pions produced in interactions of secondary nucleons is about half the size of 
that associated with the re-weighting of secondary pions. However the error 
associated with the production of the related secondary nucleons is of the 
same order. The achievable precision on the flux prediction based on thin- 
target data alone is therefore limited due to the uncertainty on the tertiary 
component of the flux. 
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Figure 3: Systematic errors of the 
[top left] and i^e [top right] fluxes at the 
far detector of T2K for the analysis de- 
scribed in Ref . |3J [7] . The beam line uncer- 
tainty combines contributions from the pro- 
ton beam, off-axis angle, target-horn align- 
ment and horn current uncertainties. The 
bottom plot shows the breakdown of the 
uncertainty on the pion multiplicity in sec- 
ondary and tertiary contributions for the i'^ 
fractional error. 



Measurements of particle emission from a full-size replica of the T2K tar- 
get have the advantage to cover at once the production of secondary particles 
exiting the target, as well as the emission of particles originating from sec- 
ondary interactions inside the target. Such measurements can be used in a 
single-step approach in which simulated yields of outgoing particles are di- 
rectly re- weighted by yields measured at the surface of the target. In this 
case, uncertainties on the flux predictions are almost entirely limited to the 
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uncertainties of the measurements. Actually, as depicted in Fig. [2} at the 
peak of the beam energy spectrum the secondary component and the tertiary 
component due to interactions in the target sum up to 90 % of the (z/g) 
flux. Hadron emission measurements with the replica of the T2K target (i.e. 
yields of charged pions and kaons exiting the target) can thus constrain up 
to 90 % of the flux prediction. 

Note that both thin-target and replica-target based approaches are nec- 
essary as discrepancies observed in a comparison of a flux prediction based 
on thin-target data to one obtained when yields of outgoing particles are 
re-weighted with the replica-target data would point to an inappropriate 
re- weighting of the secondary interactions in the target. Such comparisons 
would allow further precise tuning of the employed hadron production model. 

3. The NA61/SHINE replica-target measurements for T2K 

3.1. Experimental setup 

The NA61 detector is a large acceptance spectrometer located in the 
North Area H2 beam line of the CERN SPS. Most detector components 
were inherited from the NA49 experiment and are described in detail in 
Ref. [25]. The detector consists of a set of 5 time projection chambers (TPCs). 
Two of them, called Vertex TPCs (VTPC-1 and VTPC-2), are placed inside 
superconducting dipole magnets. The magnetic field was set to 1.14 Tm in 
order to optimize the geometrical acceptance for the T2K measurements. A 
small TPC is placed between VTPC-1 and VTPC-2 and is referred to as 
the GAP TPC. Two large TPCs, the Main TPCs (MTPC-L and MTPC-R), 
are positioned downstream of the VTPC-2, symmetrically with respect to 
the beam line. The set of TPCs is complemented by time-of-flight (ToF) 
detectors located downstream of the MTPCs. An overview of the detector is 
shown in Fig. |4] together with the definition of the NA61 coordinate system. 

The replica of the T2K target used in NA61 consists of a 90 cm (1.9 A/) 
long graphite rod of density p = 1.83 g/cm^. The downstream face of the 
target was located 52 cm upstream of the mylar entrance window of VTPC-1, 
and the target was held in position by aluminium support flanges fixed at 
its upstream end. The replica and the actual target of T2K in its complete 
environment are shown in the drawings in Fig. |5| 

A 15 kHz beam rate was used during the 2007 measurements. Due to the 
thickness of the replica target each beam proton is assumed to interact in the 
target and the trigger simply consists of selecting all beam protons by using 
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Figure 4: An example of a reconstructed p+C interaction at 30 GeV beam energy in 
the replica of the T2K target showing tracks reconstructed in the TPCs and associated 
with hits in the ToF-F detector. The incoming beam direction is along the z axis. The 
magnetic field bends the trajectory of outgoing charged particles in the x — z (horizontal) 
plane. The drift direction in the TPCs is along the y axis. 

a coincidence of various counters and vetos along the beam line (see [T3] for 
more details). In particular, the so-called SI scintillation counter provides 
timing information and triggers the data acquisition from the TPCs and ToF 
detectors. The 100 ns dead time of SI results in a 0.2 % pile up probability. 
The trajectory of each beam proton is reconstructed in a telescope of three 
beam position detectors that allows the determination of the position of the 
beam at the upstream face of the target with a precision of better than 
300 yum in both directions. 

More details on the experimental setup, detector calibration and perfor- 
mance as well as a description of the proton identification in the beam are 
given elsewhere [T3] . 

3.2. Coverage oj the T2K kinematical phase space in NA61/SHINE 

The phase space of interest for positively charged pions that exit the T2K 
target and produce neutrinos in the direction of the far detector is depicted 
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Figure 5: Technical drawing (side view) of 
the rephca target used during the NA61 data 
taking [top left] consisting of a 90 cm long 
graphite rod and aluminium support flanges. 
Drawing of the complete geometry of the 
T2K target [top right] . The overlaid red rect- 
angle represents the simplified geometry of 
the replica target. View of the T2K target 
and its cooling envelope embedded in the first 
focusing horn of the T2K beam line [bottom] . 



in Fig. |6]as a function of (p, 9), where p is the laboratory momentum of the 
pion at the surface of the target, and 9 is the angle of its direction calculated 
with respect to the beam axis. For comparison the binning used in the NA61 
data analysis is overlaid. 

The phase space of interest is divided into two kinematical regions: pions 
which exit from the side of the target with emission peaking at large angle and 
low momentum, and pions exiting from the downstream face which populate 
mainly the region of small angle and large momenta. In the T2K beam line 
the latter are not (or less) focused by the magnetic horns and are mainly the 
pions that decay to muons with momentum larger than 5 GeV/c. These 
muons are detected by the muon monitor (MUMON) located downstream of 
the beam dump and provide a spill-by-spill monitoring of the direction of the 
beam [21 ES] • The comparison of the MUMON measurements to the beam 
simulation is thus an important step in the validation of the MC model. For 
that purpose, a dedicated run was taken in 2010 in NA61 to measure precisely 
the very forward region of particle production below 20 mrad polar angle. 
In this run, the spectrometer was operated with the highest magnetic field 
configuration (of about 9 Tm), which deflected forward going particles into 
the sensitive regions of the TPCs, thus avoiding the uninstrumented region 
along the beam axis. 

The binning for the analysis in {p, 9) is driven by the acceptance of the 
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Figure 6: Kinematical phase space of positively charged pions (for 10^^ pot) exiting from 
the side of the target [left] or from the downstream face [right], and producing neutrinos 
in the direction of the far detector of T2K. The respective analysis binning of the NA61 
data with the replica of the T2K target is overlaid on top. Predictions obtained from the 
T2K beam simulation. 



NA61 apparatus. As shown in Fig. [6} it covers most of the region of interest 
for T2K. The relatively large size of the bins, ranging from 0.8 to 3.2 GeV/c 
in momentum and from 40 to 120 mrad in polar angle, is due to the low 
statistics of the 2007 data. In addition to {p, 6), data are further binned with 
respect to the longitudinal position of the outgoing particles at the surface 
of the target. As shown in Fig. |5} part of the T2K target is embedded in the 
first magnetic horn of the beam line. In this configuration and due to the 
extension of the target, the focusing properties of the horn depend on the 
longitudinal position of the outgoing particles. We investigated this effect 
with the T2K beam MC and determined that at least five longitudinal bins 
are required to obtain a prediction that does not differ significantly from a 
nominal non-binned prediction in terms of mean neutrino energy and overall 
normalisation. Five bins of 18 cm each are therefore used along the beam 
direction. An additional bin is used for the downstream face of the target. 

The acceptance of the NA61 detector in {p, 9) does not vary by more than 
10 % over the length of the target for pions exiting the side of the target. 
An identical (p, 9) binning is therefore applied to all longitudinal bins along 
the target. For pions exiting the downstream face of the target, the coverage 
extends to higher momenta. The same binning in p is maintained while a 
finer binning is used for the polar angle 9. The azimuthal acceptance of the 
detector in the x — y plane is however highly non-uniform due to the finite 



12 



extent of the TPCs along the drift direction {y axis) and the uninstrumented 
region of the detector along the beam line. This is illustrated in Fig. [7] which 
depicts the distribution of azimuthal angle (f) of the TPC tracks. For this 
reason, the NA61 replica-target data cannot be used as a direct input on 
a track-by-track basis in the T2K beam simulation for the flux predictions. 
Other suitable methods are therefore considered in Section 16.11 




3.3. Event selection and data normalisation 



As mentioned in Section [33] the NA61 beam is defined by a set of scin- 
tillation and veto counters along the beam line and the proton beam tracks 
are reconstructed in three beam position detectors. The beam tracks are 
further selected to assure that protons hit the upstream face of the target. 
The selection is based on two main cuts: the first one on the of the fit of 
the beam tracks, the second on the extrapolated position on the upstream 
face of the target. The selection rejects 32 % of the events. 

The distribution of beam particles in time with respect to the trigger time 
is shown in Fig. |8]over a 40 /is time window. Due to the relatively high beam 
intensity, about 40 % of the events include a second beam particle within 
±25 /is around the trigger time. The acquisition window of the TPCs extends 
over a maximum drift time of 50 /xs for the gas composition and drift voltages 
applied in 2007. Multiple interactions can therefore occur in the target during 
a single acquisition window. Such interactions result in so-called off-time 
tracks, i.e. tracks reconstructed in the TPCs but not associated in time with 
the beam proton that triggered the acquisition system. 
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Figure 8: Time distribution of beam particles in a 40 fis time window for single beam 
particle events [left] (~60 % of all events), and events with two beam particles [right] (^-^40 
% of all events). The beam time is centered at -300 ns. 



Since the measured yields are normalized to the number of protons on 
target, tracks reconstructed in the TPCs are associated to the triggering 
beam proton by requiring a signal in the appropriate ToF-F detector. 




Figure 9: Track multiplicity in the TPCs without (solid) and with (dashed) the ToF-F 
requirement for events with different numbers of beam particles [left] . Multiplicity distribu- 
tions normalised to the number of single beam particle events with the ToF-F requirement 
[right]. 



Actually, for the 2007 beam rate, tracks that leave a valid signal in the 
ToF-F can only have been produced in the interaction of the same beam 
proton in the target since the 100 ns acquisition window of the detector is 
much smaller than the mean distance in time between two beam particles. 
Hits associated with off-time tracks in the ToF-F detector result in overflows 
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which are cut out at the analysis leveL The effect of this cut on the track 
multiphcity in the TPCs is depicted in Fig. |9| Ahhough many beam particles 
are present in a ±25 /xs window around the beam time, the track multiplicity 
in the TPCs is consistent with that of single-interaction events once the 
ToF-F requirement is applied. 

The NA61 yields from the replica of the T2K target are thus normalised to 
the total number of protons on target which produced a valid trigger. After 
the quality cuts described above, a total of 114 885 events were selected for 
this analysis. 

3.4- Reconstruction of track parameters at the surface of the target 

Reconstruction algorithms applied for the analysis described here are 
based on those used to produce the NA61 thin-target results with the excep- 
tion that the fitting procedure at the primary interaction vertex is replaced 
by a backward extrapolation procedure to the surface of the replica target. 
The main steps of the reconstruction are: 

(i) cluster finding in the TPC raw data and calculation of the cluster 
weighted mean position and total charge, 

reconstruction of local track segments in each TPC separately, 

matching of track segments from different TPCs into global tracks. 



m 



IV 



VI 



track fitting through the magnetic field and determination of the track 
parameters at the first measured TPC cluster, 

matching of ToF-F hits with TPC tracks, 

backward extrapolation of the global tracks from their first measured 
TPC cluster to the surface of the target. 



The backward extrapolation procedure is depicted in Fig. [TOj If the ex- 
trapolated trajectory crosses the surface of the target at a certain position, 
the track parameters and associated covariance matrix are determined at this 
point. Otherwise a minimisation procedure is performed along the length of 
the target to find a point of closest approach between the track trajectory 
and the surface of the target. The track parameters are then determined 
at this point. Tracks are associated with the target if the point of closest 
approach is found within 0.6 cm from the surface of the target. This value 
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Figure 10: Sketch depicting the backward extrapolation of TPC tracks onto the surface 
of the target. The point of closest approach is determined and the track parameters p 
and 9 are calculated at this point. Only tracks for which this point lies within a distance 
of 0.6 cm around the target surface are accepted. The resolution for the different track 
parameters are also given in the figure. 



actually corresponds to the mean radial uncertainty of the extrapolation over 
the full length of the target. 

The resolution of the track parameters, p and 6, at the surface of the 
target is driven by that estimated at the first fitted TPC cluster. The latter 
strongly depends on the track topology. In order to improve the resolution, 
tracks are therefore grouped into four topologies and specific cuts on the 
minimum number of clusters on track are applied to each class. For all 
tracks a minimum number of 40 clusters is required in the MTPCs as well 
as a valid signal in the ToF-F detector. The following topologies are defined: 
the VTPC-l+VTPC-2 topology corresponds to tracks with segments in both 
VTPCs, while the VTPC-1 and VTPC-2 topologies correspond to tracks with 
a segment in one VTPC only. The GAP TPC topology corresponds to tracks 
which have measured points only in the small GAP TPC and a MTPC. 
Examples of such topologies (VTPC-2, GAP TPC and VTPC-l+VTPC-2 
from top to bottom) are shown in Fig. |4j A minimum of 40 clusters in the 
VTPC-1 is required for the VTPC-1 topology, 45 clusters for the VTPC-2 
topology, 50 clusters for the VTPC-l+VTPC-2 topology and 6 clusters for 
the GAP TPC topology. In addition, tracks are required to be reconstructed 
in a ±30 degree symmetrical wedge in the azimuthal angle with respect to 
the X-axis. The quality cuts mentioned above are used to define the detector 
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acceptance for all related MC studies in what follows. 



The resolution of p and 9 at the first TPC cluster are shown in Fig. 11 



as 

a function of momentum for the different toplogies. In particular, the GAP 
TPC tracks have their momentum measured with a maximum of 7 clusters 
in the magnetic field in the very forward region of the spectrometer. Hence 
the larger error on the polar angle and a worse momentum resolution. The 
resolution obtained after the backward extrapolation to the surface of the 
target is estimated to be = 5 cm and ere/ 9 = 6% for the longitudinal 
position and polar angle respectively. 




Figure 11: Error on the polar angle [left] and momentum resolution [right] as a function 
of momentum of the TPC tracks for data (solid) and MC (dashed). Labels refer to track 
topologies defined in the text. 

A precise knowledge of the relative alignment of the target and the beam 
is needed to reconstruct tracks at the surface of the target in bins of (p, 9, z). 
The position of the long target was first measured by surveyors. In addition, 
a procedure based on the backward extrapolation of the TPC tracks was 
developed to refine the measured position of the target with respect to the 
beam axis. For that purpose, the position of the upstream face of the target 
is used as a reference. It is actually precisely determined by the independent 
extrapolations of the TPC tracks from the downstream region, and that of 
the beam tracks from the upstream region. 
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Once the target position is known, the beam profile and radial distribu- 
tion on the upstream face are determined by extrapolating the beam tracks 
reconstructed in the beam position detectors. These distributions are shown 



in Fig. [12] together with the positions of the upstream and downstream faces 
of the target. 
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Figure 12: Profile [left] and radial distribution [right] of the beam on the upstream face 
of the replica target. The radial distribution is shown before (solid) and after (dashed) 
applying a beam track selection defined in the text below. The solid (dashed) circle shows 
the position of the upstream (downstream) face. The dotted vertical line shows the radius 
of the target. 



In 2007, the target was shifted upwards by 0.4 cm and tilted in the hor- 
izontal (vertical) plane by 5 (2.8) mrad. The hardware target alignment 
technique was improved before the 2009 data-taking period. For this data 
set, the target is well aligned along the direction of the beam (no tilt), but 
slightly shifted by 0.2 (0.1) cm in the vertical (horizontal) plane. As depicted 
in Fig. [T3| the length of the replica target is well reconstructed using the back- 
ward extrapolation procedure, which takes into account the transverse shifts 
and tilt of the target in the 2007 alignment configuration. 

In order to measure yields of outgoing particles in a configuration as 
close as possible to that of T2K (i.e. with the target aligned along the beam 
axis), beam tracks were selected to hit the target over the overlap region of 
the upstream and downstream faces, thus retaining only beam protons that 



effectively pass through the full length of the target (see Fig. 12). The effect 
of the target tilt on the yields of outgoing particles was studied over the 
analysis binning with dedicated MC simulations, and finally treated as an 
additional systematic uncertainty. 
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Figure 13: Distribution of the point of closest approach of the TPC tracks in the x—z [left] 
and y— z [right] projections after backward extrapolation to the surface of the target. The 
fact that the side of the target appears fuzzy in the vertical projection (y) is a consequence 
of the azimuthal acceptance of the detector (see Fig. [7| which is further constrained by 
the ±30 degree wedge cut defined in the text. 



The beam and target configurations in T2K and NA61 differ also by the 
beam profile on target. This could be taken into account by re- weighting 
the NA61 results with the T2K beam profile in the T2K beam MC. Due to 
the low statistics of the 2007 data such a re-weighting was not implemented. 
However, it will be applied in the analysis of the 2009 and 2010 data. For that 
purpose, the trigger hardware and software were upgraded before the 2009 
data taking. In particular, a multi-trigger acquisition system was introduced 
allowing pre-scaling of different trigger types. A certain fraction of the events 
were recorded in a configuration that defines a beam with uniform coverage of 
the upstream face of the target, and in a configuration that defines a narrow, 
centered beam. 

3.5. Particle identification 

The particle identification (PID) in NA61 relies on energy loss measure- 
ments, dE/dx, in the TPCs and the time-of-fiight that is used to compute 
the particle mass squared, m^. For each TPC track, the dE/dx is calculated 
by ordering the reconstructed clusters by increasing charge and averaging the 
distribution over the lower 50 %. For the calculation of the mass squared, 
the momentum is taken without vertex constraint and the path length of the 
track is calculated from a plane located at the center of the target along the 
beam axis to the ToF-F detector. The dE/dx and mass squared distributions 
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of the data are shown for all tracks in Fig. 14 (top panel) as a function of 
the track momentum. 
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Figure 14: Top panel: dE/dx [left] and mass squared [right] distributions for all TPC 
tracks as a function of the track momentum at the first fitted TPC cluster. Bottom panel: 
{m^, dE/dx) distributions of positively charged tracks for 40 < 9 < 100 mrad polar angle 
and 2.4 < p < 3.2 GeV/c [left], 4 < p < 4.8 GeV/c [right] momentum at the surface of 
the target. 



The dE/dx can provide an efficient PID below 1 GeV/c momentum and 
along the relativistic rise region, but is limited in the momentum region 
between 1 and 3 GeV/c where the different Bethe-Bloch curves overlap. The 
time-of-flight provides a good discrimination between pions and protons up to 
6 GeV/c. The analysis of the NA61 data with the T2K replica target is based 
on the combined PID method developed for the thin-target data analysis [12] . 
Actually, the combination of the dE/dx and time-of-flight provides a powerful 



PID over a wide momentum range. The method is illustrated in Fig. [14 
(bottom panel) which depicts how the different particles (p, vr and e) can 
be separated in the {m^, dE/dx) plane. 

A {m^, dE/dx) distribution for positively charged tracks is obtained for 
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each bin in (p, 9, z) determined at the surface of the rephca target. The data 
distributions are then fit to joint probabihty density functions (pdf) for the 
mass squared and the energy loss. Due to the independence of the dE/dx 
and w? variables, the joint pdf reduces to the product of the corresponding 
marginal distributions which are described by Gaussian distributions. The 
complete pdf is a sum of two-dimensional Gaussian distributions over four 
particle species, vr and e. For the initialisation of the fit, the resolution 

on the mass squared and the expected energy loss for each particle species are 
obtained from parametrizations of the data distributions shown in Fig. 14 as 



a function of the track momentum. The resolution on the expected energy 
loss is a function of the number of reconstructed clusters on track (~ 1/ \fN)- 
For the topology dependent cuts defined in this analysis, it is approximated 
by a constant value of 3 % due to the sufficiently large number of clusters 
on track. Independent normalisation factors are introduced for each particle 
species. Since the individual pdfs are normalised to unity, particle yields are 
given by the normalisation factors which are obtained from a two-dimensional 



log-likelihood minimisation illustrated in Fig. 15 




Figure 15: Two-dimensional fit of the data in tlie (m^, dE/dx) plane [left] and respective 
mass squared projection [right], for 40 < 6* < 100 mrad and 2.4 < p < 3.2 GeV/c at the 
surface of the target. The different components of the fit are shown with different line 
styles. 



The two-dimensional fitting procedure is applied over the full momentum 
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range of the analysis, although at high momenta when the width of the mass 
squared distribution becomes too large, the time-of-flight information cannot 
constrain the fit significantly anymore. 

Due to the low statistics of the 2007 data, a more sophisticated pdf than 
a sum of two-dimensional Gaussian distributions was not needed in this anal- 
ysis. For example, using multi-Gaussian distributions (i.e. a first Gaussian 
to describe the peak and a second one with a larger width for the tails) did 
not improve the results in terms of goodness of fit. It should also be noted 
that fits are performed in two dimensions which significantly relaxes the re- 
quirements on the pdf used to describe the data. Actually, although the one- 
dimensional Gaussian pdf's used for the dE/dx and might not describe 
tails (contaminations) of the distributions exactly, the fact that particles are 
well separated in two dimensions does not require a precise description of the 
tails in the two-dimensional case. 



4. The NA61/SHINE simulation chain 

In NA61 interactions of the incident proton beam inside the replica tar- 
get are generated as in the T2K beam simulation with the FLUKA transport 
code (the FLUKA2011.2 version was used for this analysis since the validity 
period for the FLUKA2008.3d version has already expired). The beam input 
to the standalone FLUKA simulation is based on data distributions of the 
beam divergence as a function of the position measured in the beam detec- 
tors located upstream of the target. The trajectory of each simulated beam 
track thus takes into account correlations between the position and angle of 
the beam protons. Particles exiting the target are stored and passed on as 
input to the NA61 MC detector simulation chain starting at the surface of 
the target. The GEANT3 {;28j package then propagates particles through 
the magnetic field and geometry of the detectors, and simulates physics pro- 
cesses such as particle decays. Interactions of the tracked particles in the 
detector material are simulated by the GCALOR |29] model which is also 
used for the same purpose in the T2K beam simulation. The simulated events 
are processed with the same reconstruction chain as used for the real data 
processing. 



Figure [16] shows that the employed model in the NA61 MC reasonably 
reproduces the kinematics of the tracks at the surface of the target for all the 
different topologies considered in this analysis. This is important to assure 
that the quality of the reconstruction of the track parameters is similar for 
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data and MC. Actually, the latter strongly depends on the number of clusters 
on track determined by the original kinematics at the surface of the target. 



As shown in Fig. [TT| a good agreement is obtained in terms of the resolution 
on the track parameters. 
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Figure 16: Distributions of momentum [left] and polar angle [right] of TPC tracks at the 
surface of the target for data (coloured markers) and MC (solid smoothed curves). The 
different track topologies are specified in the legend on the right plot and described in the 
text. 

As a consequence, realistic dE/dx and values are generated for the 
reconstructed MC tracks by using parametrizations of the data for the mean 
energy loss distribution and width of the m? distribution as a function of the 
track momentum (see [27] for details). 

The backward extrapolation procedure shows similar performance for MC 
and real data. An additional analysis was performed to extract yields of 
outgoing negatively charged pions in the data and the simulation. As can be 



seen in Fig. [17} good agreement is obtained between MC and data for the 
momentum distribution of negatively charged pion-like tracks after backward 
extrapolation, requirement of a point of closest approach closer than 0.6 cm 
to the surface of the target and a simple (iii^/(ia;-based PID selection to reject 
electrons. In both analyses, the efficiency of the procedure was estimated to 
be at least 98 % as a function of p, 6 and z at the surface of the target. 
The PID analysis applied to data described in the previous section is 



performed identically on the MC. Figure 18 shows the result of the log- 



likelihood fit to the simulated (m^, dE/dx) distribution in the (p, 9) bin shown 



for data in Fig. 15 
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Figure 17: Momentum distribution of negatively charged pion-like tracks [left] after 
backward extrapolation, requirement for a point of closest approach closer than 0.6 cm to 
the surface of the target and a simple dE / dx-ha,sed PID selection to reject electrons. A 
first-order polynomial fit to the ratio of data to MC [right] returns a constant of 1.046. 




Figure 18: Two-dimensional fit of the simulated data in the (to^, dE/dx) plane [left] and 
respective dE/dx projection [right], for 40 < 9 < 100 mrad and 2.4 < p < 3.2 GeV/c at 
the surface of the target. The different components of the fit are shown with different line 
styles. 
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5. Yields of positively charged pions at the surface of the replica 
target 

Yields of positively charged pions were extracted in bins of {p, 6, z) at the 
surface of the target for real and simulated data, using a log-likelihood fit (see 



Section 3.5) in the {m?,dE/dx) plane. Spectra are presented differentially 
as a function of momentum for different angular intervals, and different lon- 
gitudinal bins along the target. For simplicity, the notations duNAei/dp and 
dnpLUKA/ dp are used to refer to the data and simulated momentum spectra 
respectively, in a given angular interval and longitudinal bin. 
For data, the differential corrected spectra are defined as: 



n ri(r, ft ^\ ' 



dp Ap N^/^,,\^mp,9,z 

where Nnagi is the measured raw yield (i.e. after reconstruction and PID 
analysis) in a given angular interval and longitudinal bin for a momentum 
bin of width Ap, N^^^^ is the number of protons on target selected for 
the analysis, and the C*'s are correction factors that depend on the track 
parameters (p, 6, z). It was checked that track migration between bins is well 
below 10 % and thus the unfolding of the measured spectra is not necessary. 

Similarly, the differential spectra obtained for FLUKA with the same PID 
analysis are defined as: 

dnpLUKA _ NmC 1 TT 1 /cy\ 

Tp Ap ivf;^llc^(p,^,z) ' ^> 

where Nmc is the simulated raw yield in a given angular interval and lon- 
gitudinal bin for a momentum bin of width Ap, and N^\^ is the number of 
protons on target generated for the simulation. The Nmc raw yield contains 
part of the original FLUKA information which is reconstructed within the 
acceptance of the detector, as well as contaminations from weak decays gen- 
erated in GEANT3 and interactions in the detector material generated by the 
GCALOR model. Within the errors of the correction factors, dnpLUKA/dp is 
equivalent to the original information generated at the surface of the target 
in the standalone FLUKA simulation. 

The C* factors in Eqs. [l] and [2] include efficiencies for the reconstruc- 
tion, the backward extrapolation and the time-of-fiight detector, as well as 
corrections for the detector geometrical acceptance, pion losses (decays and 
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interactions in the detector material) and contamination from weak decays 
(feed-down). With the exception of the time-of- flight efficiency evaluated 
from the data, all the C* factors are MC based corrections. These are ap- 
plied identically to data and simulation and cancel out in the ratio of the 
data and simulated yields evaluated according to Eqs. [T] and |2} 



As will be further explained in Section 6.1, the use of the NA61 2007 
replica-target data in T2K is based on the ratio of data and simulated yields. 
Thus, only raw yields are considered in what follows. The raw spectra of 
positively charged pions are defined following Eqs. [T] and [2] as: 

dNNA61 NnA61 1 1 



dp Ap KM,e%%, 



(3) 



for the data, and: 



dNMc _ Nmc 1 
dp ~ Ap N^'c 



(4) 



for the MC. For data, the ToF-F detector efficiency, e^"^]^, is evaluated as 
a function of p and 9. Due to the ToF response not being simulated in the 
NA61 MC, e\^l^ is set to 1. Thus the time-of-flight detector efficiency is the 
only correction that does not cancel out in the ratio of real data to simulation 
and consequently it is included in the definition of the raw spectra for data. 
As an example, raw spectra measured over the most upstream, central and 
most downstream longitudinal bins, as well as the spectra measured at the 



downstream face of the target are depicted in Fig. [19] in the angular interval 
[40-100] mrad for the real and simulated data. 

Systematic uncertainties on the spectra computed via Eqs. |3] and |4] arise 
from the PID and normalisation for both real data and simulation. A sys- 
tematic uncertainty due to the time-of-flight detector efficiency is accounted 
for in the data. The systematic uncertainty associated with the PID pro- 
cedure was evaluated with the MC by comparing the pion yields obtained 
from the log-likelihood flt to the generated number of pions in the sample as 
a function of the reconstructed track momentum. The full statistics of the 
MC sample was used to estimate the uncertainty in the simulation. For the 
data, an independent MC sample with statistics equivalent to that of data 
was used. The estimated systematic uncertainty varies from 1 to 3 % for 
the MC and 1 to 5 % for the data with increasing momentum. A systematic 
uncertainty of 1.4 % was assigned to the normalisation to the number of 
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Figure 19: Spectra of outgoing positively charged pious normalised to the momentum bin 
size and number of protons on target in the angular interval [40-100] mrad for the most 
upstream [top left], central [top right] and most downstream [bottom left] longitudinal bins, 
and in the angular interval [0-40] mrad for the downstream face of the target [bottom right] . 
Error bars correspond to the sum in quadrature of statistical and systematic uncertainties. 
Smooth curves show the prediction of FLUKA2011.2 associated to tracks reconstructed 
within the acceptance of the NA61 detector (described in Section 3.4). FLUKA+GCALOR 
refers to the MC yields after reconstruction and PID analysis. 



protons on target in data. It was estimated by varying the cuts used for the 
selection of the beam tracks on target. The same uncertainty is propagated 
to the MC since the simulation of the beam tracks impinging on the target is 
based on real data distributions for the beam position and divergence. The 
systematic uncertainty associated with the ToF-F efficiency comes from the 
eventual inclusion of off-time tracks in the calculation. In order to estimate 
this uncertainty a first calculation was made using the full 50 /is drift of the 
MTPCs. Additional calculations were performed over only the first and last 
25 fis drift distances. By comparing these calculations the uncertainty on 
the time-of-fiight efficiency was estimated below 1 to 3 %. 
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The total systematic uncertainties are typically 3 to 5 %. For data how- 
ever, the overall uncertainty is dominated by the statistical uncertainty which 
is in the range of 10-15 %. 

6. Re-weighting of flux predictions with long-target data 

6.1. Re-weighting methods 

At least two different approaches based on the NA61 replica-target data 
can be followed to re-weight the predictions of the model used in the T2K 
beam MC for the simulation of hadronic interactions in the target: 

1. re-weighting factors are calculated in bins of (p, 6, z) within the T2K 
simulation. In this case weights are defined as: 

w{p, z) = Ar--,(p, e, z)/N^TAv. z) , (5) 

where N'^^^^ are the NA61 measured yields at the surface of the target 
corrected for various efficiencies, detector geometrical acceptance and 
particle losses (i.e. absolute yields), and are the yields of emitted 

particles simulated within the T2K beam MC; 

2. re-weighting factors are calculated in bins of (p, 6', z) within the NA61 
simulation. In this case weights are defined as: 

w{p, 9, z) = NZliiP, 0, z)/N^%,ip, 6, z) , (6) 

where A^^r'^gi NA61 measured yields at the surface of the target 

without any corrections (i.e. raw yields), and iV^^i are the recon- 
structed yields obtained from the NA61 simulation based on the model 
used in T2K. 

In the first approach, absolute yields are obtained by applying various 
corrections to the measured raw yields. This approach has the advantage 
that the corresponding re- weighting factors are almost model independent. 
Actually, dependencies on the model used in the NA61 MC occur only via 
several relatively small correction factors. This includes in particular losses 
due to secondary interactions in the detector material or contamination from 
weak decays that result in a maximum 5 % correction in the NA61 2007 
thin-target analysis for positively charged pions [15] . 

In the second approach, which was chosen for the analysis presented in 
this paper, there are two prerequisites: the same MC model must be used 



28 



in the T2K simulation and the NA61 analysis, and the simulated data in 
NA61 must go through the same reconstruction and PID analysis procedure 
as the real data. In this case, re-weighting factors can be calculated from 
raw yields in both data and MC since all common corrections used to obtain 
absolute yields in the first method will cancel out in the ratio. Thus we 
avoid introducing additional systematic errors on top of the large statistical 
uncertainties of the low-statistics 2007 data. However the re-weighting factors 
obtained in this way are specific to the common version of the model used in 
both simulations (i.e. if the model were to be changed in the T2K simulation, 
a new set of re-weighting factors would have to be calculated within NA61). 

Unlike thin-target based re-weighting factors which are calculated as ra- 
tios of production cross-sections, factors calculated with the replica-target 
data in both methods described above are based on yields of outgoing par- 
ticles that depend upon the beam parameters of the NA61 measurements. 
Thus, a relative re-weighting of the NA61 and T2K beam distributions is nec- 
essary when beam distributions differ significantly in the two experiments. 
Eqs. |5] and [6] should then slightly be modified to account for that additional 
degree of freedom. The final NA61 results with the replica of the T2K target 
based on the high-statistics 2009 and 2010 data sets will be obtained by us- 
ing the first approach which provides absolute particle yields per proton on 



target. As explained at the end of Section 3.4, the high statistics data will 
allow for the accounting for the relative re- weighting of the NA61 and T2K 
beams on target. 

Note that a total systematic error of typically 7 % was estimated for 
pion spectra obtained from the 2007 thin-target data [15j. Some of the 
contributions to the total systematic uncertainty (e.g. feed-down correction) 
are expected to be significantly smaller for the T2K replica-target data. Thus, 
for absolute yields of particles measured at the surface of the target, we expect 
a precision of 5 % or better. 

6.2. Application to the T2K beam simulation 

T2K beam MC predictions (based on FLUKA2011.2) can be re-weighted 
with the NA61 2007 replica-target data by calculating the re- weighting fac- 
tors defined in Eq. |6} Using Eqs. |3] and |4| these are given for each (p, 6*, z) 
bin by: 

W{p,e,z) = — , ToF ■ (7) 
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Figure 20 shows the re-weighting factors corresponding to the spectra 
depicted in Fig. 19, measured over the most upstream, central and most 
downstream longitudinal bins, as well as at the downstream face of the target. 
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Figure 20: Re- weighting factors for outgoing positively charged pions in the angular inter- 
val [40-100] mrad for the most upstream [top left], central [top right] and most downstream 
[bottom left] longitudinal bins, and in the angular interval [0-40] mrad for the downstream 
face of the target [bottom right]. Error bars correspond to the sum in quadrature of 
statistical and systematic uncertainties. 



In addition to the systematic uncertainties arising from the PID analysis, 
the normalisation and the time-of-flight detector efficiency, sources related 
to differences between the T2K target and the NA61 replica are accounted 
for in the total systematic uncertainty of the re-weighting factors. Dedicated 
FLUKA2011.2 simulations were performed to estimate the systematic un- 
certainties on the yields of outgoing charged pions due to differences in the 
replica-target geometry (i.e. contribution of the aluminium support flanges), 
alignment and density (1.83 g/cm'^ for the replica, 1.804 g/cm^ for the T2K 
target). The estimated uncertainty (within the statistical precision of the 
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simulations) was below 3 % for the differences in target geometry and density, 
while an overall 3 % uncertainty was assigned for the target misalignment. 
An additional systematic uncertainty (< 2 %) was estimated to account for 
the measured width of the beam momentum distribution which is not simu- 
lated in the NA61 MC. 

The overall systematic uncertainty on the re- weighting factors is typically 
about 6 %, with main contributions from the PID analysis at large momen- 
tum and from the target misalignment. The total error is however dominated 
by the statistical uncertainty which varies between 10 and 15 %. 



In order to use the re- weighting factors calculated with the NA61 replica- 
target data in the T2K beam simulation, a new class was implemented in 
the existing re-weighting software based on the NA61 thin-target data (de- 
scribed in Ref. |2l])- The class is implemented in such a way that either 
of two procedures can be followed to re-weight the production of positively 
charged pions: use of the thin-target data to re-weight the secondary and ter- 
tiary production in the target, or use of the replica-target data to re- weight 
outgoing pions at the surface of the target. A common re- weighting method 
is used for hadronic interactions that occur outside the target. 

For illustration of the complete re-weighting procedure, the T2K beam 
simulation was run with default beam parameters in FLUKA2011.2 and horn 
currents set to 250 kA. The prediction of the flux at the far detector re- 
weighted with the replica-target data is shown in Fig. 21 (left) together with 
the prediction re-weighted with the thin-target data. 

For the replica-target re-weighted prediction, the maximum possible er- 
rors are shown and correspond to a fully correlated 1-sigma shift of the pion 
re-weighting factors only. In the case of the thin-target re-weighted predic- 
tion, two sets of errors are shown: the first one corresponds to the total error 
which is shown in fractional form in Fig. [3} the second one corresponds only 
to the error associated with the pion multiplicity (shown in Fig. [3] as well) 
and can be compared directly to the error shown for the replica-target re- 
weighted prediction. Large uncertainties above 2 GeV neutrino energy for 
the thin-target based prediction are dominated by the error propagation of 
the kaon re-weighting. 



The ratio of the two predictions is shown in Fig. 21 (right) and indicates 



good agreement between the results of both methods. Errors on the ratio cor- 
respond to the error propagation in quadrature where only errors associated 
with the pion multiplicity are considered for the thin-target based prediction. 
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Figure 21: Re-weighted i'^ flux predictions at the far detector of T2K based on the NA61 
thin-target and rephca-target data [left] and ratio of the two predictions [right]. Details 
about the associated errors are given in the text. A linear fit to the ratio [right] is shown 
by the solid line. 

Both re- weighting methods are consistent within the uncertainties considered 
in this study for the re-weighting of the pion multiphcity. Although uncer- 
tainties are of the same order for the two approaches, it should be noted 
that in the case of the long-target based re-weighting, results were obtained 
with half the statistics of the thin-target case. The analysis of the 2009 and 
2010 long-target data will not only significantly decrease the dominant sta- 
tistical uncertainty but also some of the currently large systematics (target 
misalignment). 

The relative re- weighting of the NA61 and T2K beam distributions is not 
included at this stage of the analysis but is not expected to significantly alter 
the comparison presented here as a simple illustration of the re-weighting 
procedure. 

7. Summary and conclusions 

Precise predictions of the initial neutrino flux are needed by the T2K long- 
baseline neutrino oscillation experiment in Japan. This paper argues that 
the highest precision predictions can be reached by detailed measurements 
of hadron emission from the same target as used by T2K exposed to a proton 
beam of the same kinetic energy of 30 GeV. The corresponding data were 
recorded in 2007-2010 by the NA61/SHINE experiment at the CERN SPS 
using a replica of the T2K graphite target. 
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First, details of the experiment and data taking were described. Second, 
results from the pilot analysis of the NA61 data taken in 2007 with a replica 
of the T2K target were presented. Yields of positively charged pions were 
reconstructed at the surface of the replica target in bins of the laboratory mo- 
mentum and polar angle as a function of the longitudinal position along the 
target. Third, re- weighting factors for the model used to simulate hadronic 
interactions in the T2K target were calculated using these measurements. As 
an illustration of the complete procedure, the re-weighting factors were prop- 
agated to the neutrino flux prediction in T2K. The prediction obtained in this 
way for the z/^ flux at the far detector of T2K was finally compared to that 
obtained with a re-weighting based on the NA61 thin-target measurements. 

In the global framework of accelerator-based neutrino oscillation experi- 
ments, the paper demonstrates that high quality long-target measurements 
can be performed with the NA61 setup and that such measurements will 
lead to a significant reduction of systematic uncertainties on the neutrino 
flux predictions in long-baseline neutrino experiments. 
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